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CHAPTER  I 


INTRODUCTION 


If  a  thin  flat  plate  is  cut  from  an  orthotropic  material  (i.e., 
material  possessing  three  mntually  perpendicular  planes  of  elastic  sym- 
metry) such  that  the  plane  of  the  plate  is  parallel  to  a  plane  of  elastic 
oymnetry,  then  the  plate  is  said  to  be  orthotropic  and  it  irlll  contain  two 
perpendicular  axes  of  syrranetry.  It  is  the  purpose  of  this  dissertation  to 
consider  the  state  of  stress  in  a  large  rectangular  plate  of  this  type 
having  a  small  elliptic  disk  of  a  different  orthotropic  material  at  its 
center.  Both  the  axes  of  symmetry  of  the  plate  and  the  axes  of  elastic 
synnetry  of  the  disk  are  taken  parallel  to  the  edges  of  the  plate.  The 
dimensions  of  the  disk  are  considered  extremely  small  relative  to  those 
of  the  plate.  As  shown  in  figure  1,  the  center  of  the  disk  will  be  taken 
as  the  origin  of  the  coordinate  system  with  the  major  and  minor  axes  of 
the  elliptic  disk  coincident  with  the  coordinate  axes.  The  plate  will  be 
subjected  to  a  uniform  tension  S  along  two  opposite  edges  and  the  result- 
ing stresses  in  the  plate  and  in  the  disk  will  be  determined.  Throughout 
this  discussion  it  is  of  course  assumed  that  strains  are  small  and  wUl 
remain  within  the  limits  of  perfect  elasticity, 

A  second  objectiTe  of  this  dissertation  is  to  show  that  expressions 
for  stresses  of  several  special  problems  in  elasticity  can  be  obtained  as 
limiting  cases  of  this  more  general  problem.  Some  of  these  have  already 
been  solved.  Professor  C.B.  Smith"^  treated  the  problem  of  an  infinite 

n4.+,^Ji;  Effect  of  Elliptic  or  Circular  Holes  cp  the  Stress 

Dlstrtbutlon  in  Flates  or  >vooa  or  KLywood  Considered  as  Urthotrooie  

jv-a^enaiSt  u.  a.  i-arest  nroauota  Uhoratmy  Umpfrr*.  il'iii 
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Figure  1,  —  Elliptic  Disk  of  Orthotropic 
Material  in  an  Orthotropic  FLato 
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rectangular  orthotropic  plate  with  a  small  elliptic  or  circrxlar  hole  at 
its  center,  and  Dr,  A«  J.  Owens    treated  the  problem  of  an  infinite  ortho- 
tropic plate  with  a  small  elliptic  hole  at  the  center  filled  with  a  rigid 
elliptic  disk,  Dr,  J.  C,  Ehglish^  treated  the  problem  of  an  infinite 
rectangular  orthotropic  plate  containing  a  circular  hole  filled  by  an 
isotropic  material  and  then  showed  that  the  problems  of  Smith  and  Owens 
were  obtained  as  special  cases.  Oace  the  stresses  for  the  more  general 
problem  treated  here  are  obtained,  no  great  difficulty  is  Imposed  in 
deriving  a  few  special  cases  although,  in  some  instances,  the  forms  in 
which  the  stresses  appear  may  not  coiiwide  with  those  previously  derived 
by  other  authors.  As  an  example,  one  particular  problem  known  to  have  been 
solved  is  that  of  the  elliptic  isotropic  disk  in  an  infinite  isotropic 
plate,  L.  H.  Donnell^  determined  the  stress  distribution  for  this  case  but 
his  use  of  a  different  coordinate  system  and  method  of  attack  would  make 
a  direct  coii5)arison  with  the  results  of  this  dissertation  decidedly 
unprofitable  at  the  present  time. 

In  the  chapter  which  follows,  general  expressions  for  the  stresses 
in  the  orthotropic  plate  and  in  the  orthotropic  disk  are  derived  by  the 
use  of  ccwplex  functions.  The  procechxre  entails  a  discussion  of  the 
general  equations  of  plane  stress  for  orthotropic  materials,  application 


i^uiiomcioj-t;B,    voj.,   iA,   X5/i>X,  pp.  ^2^-33 

-^Julius  Clyde  English,  "Stress 
Orthotropic  Flates"  (unpublished  Pli.D, 
Florida,  1952). 
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of  mapping  functions,  and  finally,  derivation  of  the  stresses  from  series 
"Which  are  subjected  to  suitable  restrictions  and  boundary  conditicms. 

Chapter  HI  is  reserved  for  a  treatment  of  special  problems  irtiose 
stresses  can  be  obtained  by  specializing  the  solution  obtained  for  the 
problem  considered. 


CHAPTER  n 


SOLUTION  OF  THE  PROBLEM 
General  Equations  of  Plane  Stress 

The  solution  of  this  elasticity  problem  involves  finding  general 
expressions  for  the  stresses  such  that  the  state  of  stress  at  every  point 
in  the  object  under  consideration  may  be  specified.  These  stresses  are 
obtained  as  derivatives  of  general  stress  functions. 

Since  the  stress  components  vary  over  the  area  of  the  plate  (and 
disk),  when  we  arrive  at  the  boundary  of  the  elliptical  disk  the  stress 
components  in  the  plate  must  be  in  equilibrium  with  those  in  the  disk  so 
that  the  forces  in  the  plate  may  be  regarded  as  a  continuation  of  those 
in  the  disk. 

The  mathematical  formulation  of  the  condition  for  compatability  of 
stress  distribution  with  the  existence  of  continuous  functions   u   and  ▼ 
defining  the  deformation  results  in  a  differential  relation  called  the 
condition  of  compatability.  If  the  conditions  of  compatability  are  satis- 
fied then  the  existence  of  displacements  corresponding  to  a  given  set  of 
strain  con^jonents  is  guaranteed. 

From  the  above  discussion  it  follows  that  three  sets  of  conditions 
must  be  satisfied,  therefore,  if  the  function  is  to  give  us  the  true 
stresses,  these  conditions  being: 

1,  the  equations  of  equilibrium 

2,  the  conditions  of  compatability 
3«  the  boundaiy  conditions  , 

Now  the  stress  at  smy  point  in  a  given  material  is  completely 
determined  by  the  six  canponents  of  stress  at  that  pointj  i,e,,  three 


6 


normal  stresses  and  three  sheajring  stresses.    Similarly,  the  strain  at 
any  point  is  con^jletely  determined  by  the  six  con^jonents  of  strain;  i.e., 
three  extension  strains  and  three  shearing  strains.^  For  the  coii?)onents 
of  stress  and  strain,  Tdmoshenko's^  notation  ■vriLll  be  nsed.    The  stress 
conpcnaents  are: 

Ci*        (rP»  ^xy  '       tj^*  K 

If  u,  T  and  w  are  the  components  of  displacement  of  a  point  in  a  body 
at  a  position  iihose  coordinates  are  x,  y  and  z,  the  components  of  strain 
are  defined  as  foUomt 

r_<^v,  ^_c^Tr 

(1) 

^^-Jx  +  Jy  Tec  =^4-^ 


Assuming  that  the  coordinate  planes  coincide  trith  the  planes  of 
elastic  symmetiy,  the  components  of  stress  and  strain  for  orthotropic 
materials  are  related  by  the  following  equations 

^The  reader  nill  find  a  mmnarisation  of  most  of  the  symbols  and 
their  nwanings  as  used  throughout  this  paper  in  the  Appendix, 

^S,  Timoshenko  and  J,  N,  Goodier,  Theory  of  Elasticity  (2d  ed.). 
New  Yorkt  McGra^r  Hill  Con^jany,  1951,  

7 

'H.  W.  March,  Stress-Strain  Relations  in  Wood  and  Plywood 
Considered  as  Orthotropic  Materials,  l^adlsoo,  IVisconsini  United  States 
l-orest  mxtucts  laboratory  Heport  No.  Rl503,  19Ui,  pp.  2-3, 
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(2) 


xy  j/yz  M  zx 


■where  Ej^  ,  Py  ,       are  Young's  moduli  in  the       ,  y-  and  z-directions 
and  1^  is  Poisson's  ratio.    In  particular,  \)^  is  the  ratio  of  contracticm 
parallel  to  the  y-axLs  to  the  extension  parallel  to  the  x-axLs  associated 
Tfith  a  tension  parallel  to  the  x-axLs.        is  the  modulus  of  rigidity 
(modulus  of  elasticity  in  shear  )j  more  specifically,  y6/__  is  the  modulus 
of  rigidity  associated  nith  the  x-  ,  y-directions. 

For  problems  of  plane  stress,  equations  (2)  reduce  tot 

(3) 

i.e.,  Hooke's  Law  becomes 


8 


Tihere 


€j  =  «12  "22  '^7 


■ll  -  -i^   '  '22    -  -i-  ' 

^  ^  (5) 


Pf  oblems  involving  orthotropic  materials  can  be  solved  by  using 
a  stress  function  irtilch  involves  tiro  con^jlex  functions  each  of  -which  is 
a  function  of  a  single  complex  variable.    Professor  Smith^  has  shovm  that 
the  stress  distribution  in  an  orthotropic  plate  In  a  state  of  plane  stress 
can  be  obtained  from  a  stress  function  irtiich  is  the  real  part  of 

^    X  F3_(x4-ioc£y)  ^-  FgCx      ipdy)  (6) 


ufaere 
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€   -  (■22/"ll)*    -  (VEy) 
K    =[i(»ii»i2^(28i2  +  »33j 
=[K  +-(K2-l)2ji 

/(?  zg    -  (K2-l)2]i- 

It  will  be  convenient  to  make  a  slight  change  in  the  form  of  equation  (6) 
and  take  the  stress  function  as  the  real  part  of 

^  -   F3^(x-»-io<^y)      FgCx -f  iOfgy)  (8) 

g 

ra  X           ^'  Smith,  "Effect  of  Hyperbolic  Notches  on  the  Stress 
I^stribution  in  a  Wood  Plate,"  Quarterly  of  Applied  Mathematics.  Vol.  VI, 
19k9,  pp.  Ii52-U$6,   » 
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or  as  the  real  part  of 
■where 

The  quantities  and    0(2  are  real  qtiantities  that  depend  upon  the 

elastic  constants  of  the  plate. 

The  stress  components  are  then  given  by  the  real  parts  of  the 
expressions 

(TV  =:      ^Fi(x+ioLy)  -f-   il.  F  (x i  ocy)  (10) 

In  accordance  Trith  this  system  the  following  notation  irill  be  used 
for  the  plate  (i.e.,  exterior  to  the  disk): 

^   -  R  ^"fj^Cx  -f  i o^^y)  +  FgCx  +  iQfgy)^ 

r  1  (11) 

=  ^}^l(^)  +  FgCsg)  J 

Tfhere  j  stands  for  the  real  part  of  the  expression  in  the  bracket 

and  i 
0^1=  >      0(1=  >       €  =(s22/sii)* 

K  =:  |(3ii^22^^^2^2  "**  ^33)}  (12) 
CX  r {k  +  (k2  -  l)^"j^     ,    ^        -  (k2  -  1)2^2 
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Hooke's  Lew  being 


(13) 


-  Jx  ^  Ty  ^ 

■where  the  elastic  constants  of  the  plate  only  are  used. 

Similarly,  for  the  orthotropic  disk  "within  the  plate  the  following 
notation  Trill  be  used: 


(lU) 


I  -  r[^i(x +ia^y)  +  ?2(=^-*-ia2y)] 

■where 

ii.[j(8uy*(28^  + ijj)]  (15) 
«  =  [li  +  (F-i)*J*  ,   ,3 -[«  -  (R^- 1)*]* 

Hooke's  laif  being 

where  the  elastic  constants  used  here  refer  only  to  the  orthotropic 
material  of  the  disk. 

In  deriving  (6)  the  equations  of  equilibrium  and  the  ccmpatability 
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equation  have  been  satisfied   so  that  the  remainder  of  this  analysis  irill 
be  devoted  largely  to  the  boundary  conditions* 

For  the  problem  under  consideration  stresses  derivable  from  a 
stress  function  of  the  form  (U)  irtiich  satisfy  the  boundary  conditions 
at  infinity  and  on  the  boundary  of  the  disk  wiH.  be  sought  j  additionally, 
stresses  derivable  from  a  stress  function  of  the  form  (lU)  irtiich  satisfy 
the  boundary  conditions  along  the  periphery  of  the  disk  are  to  be  determined. 


) 


^aidj.  p.  k$3. 
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Boundary  Conditions 

The  equations  of  the  ellipse  in  figure  1  are  given  by 

/X.  =  a  © 

^  -  J/-         ©  • 
If  the  positive  direction  of  the  normal  ^  be  oriented  as  shown  in 
figure  2,  then  X' 
represent  the  x-  and  y-con^jonents 
of  the  force  per  unit  thickness  of 
the  plate  acting  on  the  element  of 
arc  cLa.  of  the  boundary  from  the 
positive  side  indicated  by  the 
direction  of  positive    1/  , 


(17) 


Figure  2.  —  Orientation  of 
Positive  Nomal  and  Tangent  on  the 
Boundary  of  the  Ellipse, 


Then  from  considerations  of  equilibidum 


(18) 


where 


and 


(19) 


One  observes  that  if  ^   represents  the  stress  function  for  the  plate. 


X3 


then 


X 


(20) 


Figure  3»  —  Resolution  of  Stress  CoB^Donents, 

It  is  clear  that  the  displacements  of  the  disk  and  plate  mast  be 
equal  at  the  boundary  of  the  disk,  and  it  will  be  shown  later  that  x-  and 
y-ccmponents  of  the  force  on  the  boundary  per  tmit  thickness  must  also  be 
equal  for  the  plate  and  disk.    Consequently,  along  the  periphery  of  the 
disk  it  is  required  that 


u  =:  fi  , 


Ay    Ay, 


Y    -  ? 


(a) 


and  for  the  plate  (exterior  to  the  disk)  the  boundary  conditions  are 
specified  to  be 


(22) 


where  S  is  the  uniform  tension  applied  along  two  opposite  edges  of  the 
plate. 
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Preliminary  Discussion  on  Mapping  Procedure 

Mapping  functions  irhich  are  analytic  functions  of  complex  ▼ariables 
will  be  used  to  solve  the  problem  of  the  orthotropic  plate.    In  order  to 
show  the  general  procedure  involved,  a  brief  digression  will  be  made  from 
the  original  problem  in  considering  a  sia^jle  mapping  of  an  ellipse  in  the 
a-plane  to  a  unit  circle  in  the  ^  -plane* 
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z-plane  ^  -plane 

Figure      —  Ellipse  Mapped  into  a  Unit  Circle. 

As  is  customary  in  complex  tnalysis,  set  «  =  x  +  iy  where 
On  the      -plane  this  ellipse  is  mapped  Into  the  unit  circle 


i«e.. 

Hence 


^  =  t6  where   r  =c  1  j 
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so  that 

z  ■=  X  -+•  iy 
Solving  for  ^  in  terms  of  z  it  follows  that 

i  t  ii^- (c-"- j-^) 


which  is  a  double-valued  function  of  z  and  hence  contains  tiro  branches, 
The  two  branch  points  are  found  by  setting  the  radical  equal  to  zerO) 
which  gives 


2=  ±{IZI^  . 


These  points  are  readily  seen  to  be  the  two  foci  of  the  ellipse,  and  it 
is  possible  to  take  the  line  on  the  x-axis  joining  these  two  points  as 
the  branch  cut# 

On  the  x-.axis,    y  jz  0,  so  that 


Hence 


and 
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z-plane  ^  -plane 


Figure  5«  —  Branch  Cut  of  Ellipse  Mapped  into  a  Circle. 

As  a  point  moves  along  the  x-axis  between  the  branch  points 
of  the  ellipse  in  the  z-plane,  a  corresponding  point  moves 
around  a  circle  with  center  at  origin  inside  the  unit  circle  in  the 
^  -plane;  the  radius  of  this  circle  being 

To  insure  that  ^  is  a  single-valued  function  of  %  (i.e.,  to  map  ^  Into 
the  region  outside  of  this  circle)  it  irill  be  necessary  to  have 


It  should  also  be  noted  that  as  a  point  moves  around  the  ellipse  in  the 
z-plane,  a  corresponding  point  moves  around  a  unit  circle  in  the  -plane. 
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Mapping  Applied  to  the  Qrthotropic  Disk 

By  a  procedure  similar  to  that  in  the  preceding  section,  a  mapping 
of  the  boundary  of  the  orthotropic  disk  is  effected. 
In  the  Z]^-plaiie  let 


o  .  •  5 

2,  =  oc  +  co^€/^ 


-'X+l'CX,^  (23) 

•where  0.  ,  4,=  oi^Jt-^^Q  is  the  bormdary  of  the  disk.    To  map 

a 

this  curve  into  a  unit  circle  on  the  ^^  -plane  let 


^,  =  e  ^  =  1  -h  t  n  . 


Then 


and 


or 


4 


where 


(2U) 


^.r:^-c^-^  (25) 
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Hence 


and 


or 


1  irtj^^C^^^ 


The  function      has  branch  points    z^^s  i/a^  -Scfb^    •    Hence  the 
branch  cut  will  lie  on  the   k^^-axis  if  o(-^  <  a  »    and  on  the  ^--j^-axis 
if    S<3^b  >  a.    If  Sf'^b  =s  a,  the  mapping  problem  is  trivial. 

Consider  the  case  ^^b^  a*    It  is  seen  that   f^-==.0    on  the 
branch  cut  so  that 

from  -which 


and 


Hence  the  branch  cut  in  the  ZQ^-plane  maps  into  a  cirele  whose  radius  is 
less  than  that  of  the  unit  circle  In  the    '^i--p'ljax>B,    To  map  into 
the  region  outside  of  this  circle  it  will  always  be  necessary  to  have 


19 


s 


la.'-i.-t'-  / 


• 

H— 

J  ' 
^  / 

• 

€  -plane 

Figure  6.  —  Mapping  of  Orthotropic  Disk  into  an  Annular  Region 

For  the  case  Q(-^b  >  a  »  the  branch  cut  trill  lie  on  the  yj^-axiB 
•where  x-j^z:  0,    5y  following  the  same  procedure  it  is  readily  seen  that 
the  ellipse       =:      +  iy^   maps  into  the  unit  circle       =  e"^^  and  the 
branch  cut  in  the  Z]^-plane  m^s  into  a  circle  in  the    ^^-plane  ndth  radius 

Then  to  ensure  mapping  "the  outer  part  of  this  circle  it  -will  be 

required  that 


Similar  results  ■will  be  obtained  by  carrying  through  the  same 


arguments  for  the    S^-plane  in  irtiich 


•  e 


(26) 


T^re         =:  a  cos  0  ,       y^  =  oj^gb  sin  Q  . 
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Then 


or 


■vrtiere 


(27) 


(28) 


The  ellipse  in  the  *z^-^lsjaB  then  transforms  into  a  unit  circle  in  the 
^-plane  with  the  final  results  being 


X 


■where 


depending  on  -whether  oCj^  <  a   or    <^b  ;>  a  , 

Hence  to  ensure  suitable  mappings  of  the  regirais  in  question  it 
■will  be  necessary  to  use 


and 


(29) 


where  it  shall  be  required  that  the  relations 


always  hold. 


(30) 
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Continuity  Conditions 

In  solving  ttro-dimensional  problems  for  orthotropic  materials, 
the  stress  function  (6)  T»hich  is  a  function  of  two  complex  functions  is 
obtained  by  the  aid  of  two  mapping  functions j  for  the  elliptic  disk  these 
mapping  functions  are  y  ^   and     ^2  • 

It  has  been  shovm  that  the  boundary  of  the  ellipse  in  each  of  the 
8^-  and  S^-Planes  can  be  mapped  into  tmlt  circles  in  the     f^-  and 
planes  respectively}  and  furthermore  the  branch  cuts  in  the    8^-  and 
Sg-Planes  can  be  mapped  into  circles  irLth  radii  less  than  the  unit  circle 
on  the  and   f^-^plaxieo .    The  elliptical  regions  in  the  8^-   and  the 

Sg-planes  are  then  mapped  into  annular  regions  on  the  and  ^2-Plane8, 

The  ring-shaped  regions  into  irtiich  the  ellipses  and  their  branch 
cuts  are  mapped  by  the  functions     ^   and    5     suggest  the  use  of 
Laurent  series  in   J-^  and  purpose.    Accordingly,  it  will  be 

convenient  to  take 

T»iiere   ^   and         are  con^^lex  coefficients  to  be  determined  by  the 
boundary  conditions. 

Now  by  virtue  of  equations  (10)  and  (31)  the  stress  components 
for  the  disk  can  be  written  as 
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fjvjjl.'f.     2'"'>  X*"^ 


/ 

(32) 


To  etisure  continuity  of  the  stress  function,  a  relationship  between 

the  coefficients   Xi    and  1      ,  and   8     and   S  „   must  be  obtained. 

"  -n  n  -*i 

If  it  is  assumed  that  <5'-b<a  ,  the  branch  cut  will  lie  along  the 
Xj^-axis,  and  in  the   ^^lane  we  -will  have 


I    r-  and         —  «  . 


so  that 


23 


As      approaches  the  branch  cut  from  above. 


and  as  %^  approaches  the  branch  cut  from  below 
If  is  continuous,  then  it  is  required  that 


Hence 

/4  A.  ^  /I,  A 

^  A", 

Sliiiilarl7,  if  <  a   then  on  the   '^^-plsne  one  dbtains 


(3U) 


irhere  ,  ^^^^  ,  ,  are  constants  for  the  disk  analogous  to  (2^) 
and  (28), 

Now  if  4fib  >  a  ,  the  branch  cut  ifill  lie  along  the  y^-cudsj  then 
as      approaches  the  branch  cut  from  the  right  side 

where    -%  ^0=^%. 


sr 


Aa  8^  approaches  the  branch  cut  from  the  left  side 


.6© 


If  is  continuous  as  '^■^  approaches  the  bi^ch  cut  from  the  right 
and  from  the  left,  then 


-at 

so  that 


-  4.  '^J  or    d„  =  4 

And  similarly  for  the    ^  ^-'V^-^^ 

these  results  being  identical  with  (33)  and  (3li)« 
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General  Equations  of  Displacement 
For  the  plate  (exterior  of  disk) 

from  nhich  it  irill  be  convenient  to  set^^ 


where,  from  equation  (.29) »  it  is  seen  that 
and 

Substituting  the  derivatives 


and  -yU,  (36) 


one  nor  obtains 


Whereas  it  Tras  found  convenient  to  represent  the  interior  of  the 
disk  by  a  laurent  series,  it  will  now  be  desirable  to  represent  the  strese 
functions  exterior  to  the  disk  by  series  of  negative  powers  in  ^  inasmuch 
as  the  stresses  in  the  plate  at  00  must  approach  constant  values. 
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"      Jx^        '^l    Jj,'    ^    Jit  J 

Taking  the  partial  derivative  of  ^  -with  respect  to  y. 


Now  from  (1)  and  (U) 
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Therefore 


and 


(U2) 


Integration  of  (U-)  and  (U2)  yields 


Now  from  (1),  (U)  and  (10) 
But  from  (li3) 


is 


Inhere   u^(y)    and   v.(x)   mean     ^f^and  respect ive]y»  so  that 


i.e.. 


in 


The  first  bracket  of  (i;6)  reduces  to 
Multiplying  by  ^/ this  becomes 
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T»here  it  is  recalled  tram.  (12)  that 
similarly,  the  second  bracket  in  (U6),  Tia.t 


simplifies  to 

Replacing  the  brackets  of         with  these  expressions  and  then  equating 
(Ui)  and  (U5)  reveals  that 

i*e,> 

since  the  left  hand  member  in  this  equation  is  a  functlcm  of  7  <»ly, 

while  the  right  hand  member  is  a  function  of  x  only,  tib  are  lead  to  conclude 

that  both  meinbers  must  be  equal  to  some  constant,  call  it  Cj  then 


and  ^ 


1r 
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Hence 


and 


Integrating  (U?)  and  (U8), 


AT 

that  is 


Similarly  for  the  interior  of  the  disk  ure  obtain 


(50) 


The  linear  functions  of  x  and  y  in  equations  (U9)  and  (^0)  can  be  disposed 
of  by  fixing  the  origin  and  eliminating  rotation  and  translation  with 
respect  to  that  origin.    If  this  is  done,  the  eapressions  for   u  ,  t  ,  fl 
and        reduce  to  the  terms  within  the  respective  braces j  i.e.. 
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Determination  of  Goefficienta 

The  complex  coefficients  in  the  series  representing  the  stress 
functions  can  now  be  determined  by  inqsosing  givmi  boundary  restrictions 
on  the  disk  and  on  the  plate. 

Since  a  uniform  tension  S  is  acting  upon  two  opposite  edges  of  the 
infinite  plate,  and  since        exhibits  the  sane  behavior  at  infinity  as 
s,  then  it  is  seen  from  (22),  (37),  (39)  and  (l^O)  that 


and 


Noir  from  (52) 


A,  A 


X 


and  from  (53) 


(52) 


(53) 


i'-^'^      ^      (A.  '^•^  / 


A, 

TThere  the  primes  denote  the  real  parts  of  the  con5)lex  coefficients. 
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Therefore 


and 


6:  = 


A, 
5  A. 


From 


D'  -  Aj,  


A|  A;t 


(55) 


Therefore 


//     -a,  Ax_  -j'^ 

5i  —  ^'  (56) 

nhere  the  double  primes  denote  the  imaginary  parts  of  the  conqplex 
coefficients* 

Now  at  the  boundary  of  the  disk  the  ellipse  has  been  mapped  into 
a  unit  circle  so  that  here 

y,=  f^  =  f.  =  i  =  e  ^  (57) 

Making  this  substitution  in  the  first  condition  of  (21),  i.e.,  setting 
fl      u  yields 

ft[C-  4,-* J,,^- /I. e  +  (-«.^.,+^:t)^  &,  e  j 


where  the  general  expressions  for  u  and  €  were  given  in  equations  (51), 


The  preceding  equality  can  be  written  as 


where 


C,  =  C-  c^.*  -i-"    ^'x)      ^      4  =  C-  +  A) 

Similarly  from  "9  =  v   is  obtained  the  relationship 

•where 

Dividing  out  a  two,  equations  (58)  and  (6o)  can  be  Tiritt«n  as 


and 
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frcBi  irhich  it  is  seen  that 

•where   n  =  0,  1,  2,  .  .  ,     and        =       =  0      for     n  ^  2  . 

Now  since  ^(^)^     and         -   \{^)^     by  equations  (33) 

and  (3U),  then  (6U)  and  (65)  become 


for  the  Talues   n  -  0,  1  , 

For  n  ^  2,  these  equations  reduce  to 


and 


(66) 


(67) 


(68) 


=  +  ^  ^4  ?»,  (69) 
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Now  on  the  boundary  of  the  disk  the  real  part  of  the  function  ^ 
satisfies  the  relations 


T»here   ds    denotes  an  element  of  the  boundary  and  £y  ds     and   ?p  ds 
represent  the  x-  and  y-con^jonents  of  the  force  per  unit  thickness  of  the 
disk  acting  on  the  element    ds   from  the  positive  side  as  indicated  by 
the  direction  of      ,  the  noiroal  to  the  boundary  (cf,  figure  3). 
Similarly,  in  the  plate 

But  the  boundary  of  the  ellipse  has  been  mapped  into  a  unit  circle  so 
that 

or 

±fil\± -y       and -sL/'^U  ^ - y' 

Thus  to  preserve  continuity  of  the  function  In  going  from  the  plate  to 
the  disk  the  following  relationships  are  required  at  the  boundaryi 


cJieK  <^7a  J    oi&\^'*-  I  (70) 

and 


(71) 
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Nov 


tie 


Also 


Similarly 


and 

Setting  (73)  equal  to  (72)  yields  the  equations 

[i-^K^Tyi^^  i{iT]=  forn^l  (76) 

and 

[4^>lftrj+[4,+  4(^/'J^  C  forn^2.  (77) 

Similarly,  from  (75)  and  (7li)  to  have 

-  °''k-ij+  (78) 

for  n  =^  1»  and 
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(79) 


for    n  ^  2  •     Ncfw  after  taking  cojijiigates  there  remain,  for  the  case 
n  >2,    the  four  equations  (68),  (6?),  (77)  and  (79)  in  the  unknoims  , 
^  •         ^Hi*  ^  similarly  for  8^'' »         and        j  i.e.,  for  n  >  2^ 


(80) 


and  four  similar  equations  involving  the  imaginary  coefficients  l!\ 
,        and        ,  viz., 

d+i")  K  ^   O+O  ^"  =   A"  +  8"  (81) 
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Let  the  preceding  equations  (81)  be  uritten  as 

"Where 

L=  l^i^f""  P  =  1- 

M  =  1^(82)^         Q  =1-  (82)"°  ; 

Then  it  is  seen  from  the  above  set  of  four  homogeneous  equations  in  four 
unknoims  that 

*n  =  ^„  =  A_  =         =  0        for     n  ^  2 
ix       n  — n  -n 

and  similarly  from  (80) 

C=6rr=C  =  ^^=°        for  n>:2. 
Since  the  coefficients        ,        >  etc,  for  the  case   n  =:  0   do  not  enter 
into  the  expressions  for  the  stresses  or  strains,  this  case  can  be  ignored 
so  that  there  now  remain  equations  (66),  (67),  (76)  and  (78)  for 
determining  coefficients  in  uhich  n   takes  the  value  unityj  i.e., 

[a. f  imj ^  [i^iay]  =  (A + 1, )  ^  ( 8,  ^  I) 
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Taking  the  conjugates  of  (82)  and  successively  adding  and  subtracting 
yields  the  two  sets  of  equations 


(B3) 


and 


(8U) 


Equations  (8ii)  and  (^)  involve  five  equations  but  six  unknown 
coefficients.    Therefore  ire  shall  try  setting  these  Ijaagluary  parts  equal 
to  aero;  i.e.,  set 

If  this  is  done,  then  the  stresses  irtll  involve  only  real  coefficients 
T«fhich  can  be  obtained  from  (83), 

To  solve  (83)  rewrite  it  as  followBJ 


la 


Tfhere  now 


L  =  1+  =  2aA*     ,    P  =1 

The  determinant  ^  of  the  coefficients  is 


2  ^iVa", 


A  - 


a.  p 

L 

n 

-1 

-1 

Erpanding  and  then  s\ibstituting  the  above  values  for  L,  M,  P  and  Q  leadfl 
to  the  expression 

+ a.      5!.  («,-  ■=<.)-  c,4  Cc<,  -<4)+  ai,  (c,«^-  c^cx, ) 


-I 


U2 


Now  for  the  product  AX    ky   -n^  have  the  determinant 

(A:~^d',)      M    -I  -1 


(X, 


Expanding  in  terms  of         and         and  then  making  subatitutiona  from 
equations  (5^)  leads  to  the  eocpression 

A*.  k«u  -^ijK  ^ 

Sdjnilarl7,  solving  (85)  for  ^Sj',  ^  A^^  ,  and  A  B^^   in  tems  of 

and   Ej_   irith  further  substitutions  frm.         yields  the  following  results* 

-f  (V.^ /i, (c. c/, 5,  - C,c^,        CjS. « .  -  C. 3.  ) 

+  kJrX^(cA<>^.  -^.<i'^'  -  C^y.  5:^)  j  (88) 


U3 


and 


(89) 


(90) 


Hetice  frcm  (86),  (87),  (88),  (89)  and  (90)  the  folloirLng  eocpressions  are 
obtained  for  the  coefficients. 


hh 


5  a\  r 


and 


•where 


(9U) 


(den)  =|^a[CcV4)C<^<^.^c/.c..|+a<5^.^[4  C-^Xa-'^.) 

Since  all  imaginary  coefficients  have  been  discarded,  the  express- 
ions for  the  plate  stresses  now  become  (as  seen  from  equations  39,  37 
and  UO) 
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and  from  (32)  the  disk  stresses  reduce  to 


Tfhere         and  ^2   are  given  by  (29 )  and  and  are  simllap 

erpressions  for  the  disk. 
Nor  we  hav* 

b7  equations  (33)  and  (3li)j  and  if  ^   and         are  defined  as 
then  ulth  these  substittitions  the  plate  and  disk  stresses  reduce  to 
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(98) 


(99) 


Trtiere        ,  S^'  ,       ,       ,  A^j^  and  B^^^   are  defined  by  (91),  (92), 
($S),  (93)  and  (9U)  respectively. 

Equations  (98)  and  (99)  represent  the  complete  solution  of  ths 
plane  stresses  for  the  problem  of  an  elliptical  disk  of  orthotropic 
material  intoedded  in  an  infinite  rectangular  orthotr<^ic  plate  which  ia 
subjected  to  a  tension  (or  con^sression)  S  on  two  opposing  edges. 
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It  Tiill  be  achrantageous  in  the  following  sections  to  have  the  disk 
stresses  (99)  refwritten  irith  the  stibstittitions  (91)  and  (92).  These 
results  folloirt 

5a  J- (a.  C.o^,- A,  qo^,  +A^C^oa  -  C,  cX^  ) 
a  (A,C,o/.(y^-A,C,o/,o^,-»-AjC^o/,ot,->4^>^/^u) 


(100) 


-  <?,C^    (  A,  -  A.o/j^i    +  A^<W;iO^,  -Aao/,^^) 

-  J-Ca.  CjO',  -  A,  c^ot,  -h  Aj^c,      -  Aj,c^o<^) 

It  is  noted  that  the  shear  stress  becomes  zero  since  the  quantities   o(,  , 
oij^  ,         and   Ej^  are  all  real. 


CHAPTER  in 


DETERMINATION  OF  STRESSES  FCR  SPECIAL  CASES 


Circular  Qrthotropic  dsk  in  an  Qrthotroplc  Plate 


It  will  now  be  of  interest  to  reduce  this  solution  to  several 
special  cases.    Consider  first  a  circular  disk  of  orthotropic  material 
instead  of  the  elliptic  disk.    In  this  case  it  is  merely  necessary  to  set 
b  =:  a  f  which  results  in  slight  sinpllfication  of  the  preceding  express- 
loos*    In  particular 


(102) 


/ 

^.1  = 


5cL 


U9 
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and  -  ^ 


(den/r::|^[(c.-C.)C,;/,o<._c^,o<,^ 

+   '^.jcc.-- c.)(S,  A.  -  3.  5,)  J  j  (105) 

Bquationa  (102),  (103),  (lOU)  and  (105)  are  no«r  substituted  into  tha 
stresses  (98)  irtiich  are  further  modified  by  setting 

Aj,=  ^0  +  <^u)  X,z;aC»-^-»)  (106) 


$1 


and 


^2 


^2  -  a2(l-o<i)  . 


(107) 


Stresses  for  the  circular  carbhofcropic  disk  reduce  to 


(108) 
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Isotropic  ELliprbic  Disk  in  an  Qrthotropie  Plate 

For  the  case  in  nhich  the  disk  material  becomes  isotropic,  it  is 
necessary  that 

o 

since  the  elastic  properties  are  independent  of  direction.  Hence  £"  =  1 
so  that  <X,  =  o(     and  c^j  =  /3  . 

e  J-  '  ^  11  ^ 

Nor   o(»  and  /(?  as  used  by  Professor  Smith     have  the  values  6 
and  e    respectively.    The  function 

Twas  obtained  from  the  compatability  equation 

T*here        =67  and  ,    lCp   -^CV  -'•C 

are  roots  of  the  auxiliary  quadratic  equation 

in  ifhich 

and        ^  =r  —  €  or      —  Q, 
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Consequently  the  folloirijag  relationships  exist  for  the  isotropic  elliptic 
disk  in  nhich  <^  tends  to  zero. 


^  o 

E 


o 


5U 


0        2-*/  / 
For  the  relationships  above        ,    ^  *    •'^-1  »    ^   ^-1  ^® 

indeterminant  form   0/0  •    Evaluating  by  differentiating  both  nonerator 
and  denominator  -with  respect  to        ,  and  then  letting    (p  z:.  0  and 
€  =1  ,  it  follows  that 


+-^[Cxrt-A«.)Cc.-C,3+5«-[(A.-A,X^,"<.-<4«.) 


(110) 


Si 


irtiere 


(111) 


55 


and 

1, 


 JeU  <^J^\cJ^ -  clJ^K ^^^^ 


(113) 


- 


(nil) 


Ejcpressions  (55),  (113)  and  (Hit)  substitnted  into  (98)  plus  the  stresses 
(110)  and  (111)  nor  give  the  cc3n5>lete  solution  for  the  case  of  an 
elliptic  isotropic  disk  in  an  orthotropic  plate. 

It  "vrill  be  of  further  interest  now  to  reduce  this  solution  to  the 
case  of  the  circular  isotropic  disk  in  an  infinite  orthotropic  plate. 


circular  Disk  of  Isotropic  Material  In  an  Qrthotroplc  Plate 

For  the  special  case  of  a  circular  Isotropic  disk  in  an  orthotropic 
plate  we  follow  a  previous  procedure  by  setting  b  r=  a  «  -where  it  should 
be  noted  that  both  A,  and  Aj^  contain   a   and  b« 

Rewriting   (T^   and  from  (110)  and  (111)  with  the  substitutions 

=  a  +  0(-]h  and  A2  =  a  +  O^gb  and  then  letting  b  a  yields 
the  following  results t 


where 


(115) 


(117) 
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Also 


(118) 


and 


(119) 


For  the  circular  Isotropic  disk  in  an  orthotropic  plate,  the  constants 
and    Bj^   take  the  valTies  (102),  and  equations  (98)  are  modified  by 
the  relationships  (106)  and  (I07), 
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Elliptic  Hole  in  an  Qrthotropic  Plate 

If  the  orthotropic  plate  contains  an  elliptic  hole  instead  of  a 
disk,  then  all  the  disk  constants  reduce  to  zero,  riz.f 

S=  J/  -  0 

This  substitution  makes     ^    —   0%-  —  ^Tw  —  0    as  seen  frcjm  equations 
(110),  (111)  and  (112).    Then  fi-om  (55)»  (112),  (U3)  and  (llU)  the  plate 
constants  are  found  to  reduce  to 


These  constants  substituted  into  (98)  now  provide  the  stresses  for  thia 
particular  situation. 

A  comparison  of  these  results  Trith  those  obtained  by  Smith'''^  may 
not  at  once  appear  identical,    Hovrenrer  a  slightly  different  mode  of 
attack  may  clarify  the  situation. 
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By  the  restilts  of  Chapter  II,  eqviation  (3^)  becomes 
SubstitutiJig  values  tvm.  (120)  now  yields 

M  _  /?f  _5[fa..,i);!.i^  t  ^i^^ 


Noir  since 


2zi  =  (aH-O^Lb);^  +  (a  -  g<ib) 
282  =      -^0^^)^  +  (a  -  0^2^^ 


then  equation  (121)  further  reduces  to 


and 


^  of  SA 


(la) 


(122) 
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Similarly, 


5^  [i^.^g-f-*^.^) 


A  con^iarison  of  equations  (122)  and  (123)  "with  the  results  of  Smith 
ncnr  shows  them  to  be  the  same  so  that  the  subsequent  differentiation 
BBist  lead  to  his  stresses  for  the  orthotropic  plate  nith  an  elliptic 
hole* 


(123) 
13 
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Circular  Hole  in  Orbhotropic  Plate 


If  the  hole  is  circular  instead  of  elliptic,  equations  (120)  are 
modified  by  setting  b  ss  a    so  that 


1     «A  C^/- 


Then  equations  (98)  for  the  orthotropic  plate  iiith  a  circular  hole  can 
be  written  as  / 


irhere     AL=a(l  +  0^)    ,       A2  =  a(l  +  O^g) 


and  W. 


'1   and         are  given  by  (107). 
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/      /     /  / 

Further  substitutions  for        ,  B^^  ,  A^^    and   B_-j^   yield  the  expressions 


(I2li) 


As  -was  seen  for  the  elliptical  hole,  more  convincing  results  for 
the  circular  hole  can  be  obtained  by  a  development  from  and  , 

l.e..  ,  ,        ^  "^t 
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iriilch  reduces  to 


Similarly, 

a 


^  ft 


6a^ 


(125) 


(126) 


Apporopriate  differentiation  of  (12^)  and  (126)  then  lead  to  the  same 
•tresses  as  dertred  b7  Professor  Sndth"'^, 

As  a  further  check,  consider  the  stresses  at  the  point  , 

t 

•g  ~  i  0(2*     ^  irtiich  case 


\  -  ia 


Wg  =  ia 


-  ia(l+  C<2) 


Substittttion  of  these  values  into  (12ii)  lead  to  the  expressions 

n:^^=  o. 


lii 


Ibid.,  p.  9, 
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If  the  plate  itself  Is  Isotropic  such  that 

CKi  -  0(2  =  ^ 
then  we  have  the  irell  known  resxilts"^^ 


<7J  :=  3S 


Similarly,  for  the  point   Z]^  -  a  ,    Z2  —  a     it  follows  that 

-H         =  a(l  +  Uj^) 

and  the  stresses  (12li)  become 


-  0  • 


For  the  isotropic  plate  these  reduce  to 


„.,,  .  ^^'^  tee:  C,  T.  Wang,  Applied  EUatlelty.Mew  York;  McOrair 
Hill  Conpany,  1953,  p.  62.   ^ — " 
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Svmmary 

In  brief  sMnmarization,  this  dissertation  provides  expressions 
for  stresses  in  the  following  elasticity  problems; 
1,   An  elliptic  disk  of  orthotropic  material  in  an  orthotropic  plate. 
2»    A  circular  disk  of  orthotropic  material  in  an  orthotropic  plate* 
3.    An  elliptic  disk  of  isotropic  material  in  an  orthotropic  plate, 
li.    A  circular  disk  of  isotropic  material  in  an  orthotropie  plate, 

5,  An  elliptic  hole  in  an  orthotropic  plate, 

6,  A  ciixrolar  hole  in  an  orthotropic  plate. 

For  the  reader  Tiho  is  inclined  to  investigate  farther,  four  other 
cases  could  be  derived,  yiz,» 

7,  A  rigid  elliptic  disk  in  an  orthotropic  plate. 

8,  A  rigid  circular  disk  in  an  orthotropic  plate, 

9,  A  rigid  elliptic  disk  in  an  isotropic  plate. 
10,   A  rigid  circular  disk  in  an  isotropic  plate. 

In  the  preceding  cases  of  the  rigid  disks  it  is  important  to  note  that 
Toung'8  moduli  for  the  disk  materials  approach  Oo  ijhich  irould  in^^ose 
other  obvious  modificationa  on  the  disk  constants. 

As  further  exercises,  the  reader  may  ulsh  to  reduce  iteano  1  -  U 
to  the  situations  in  nhich  the  plate  becomes  isotropic. 


APPENDIX 
General  Notation 

a  f  b  lengths  of  send'«iajor  and  seml'^ilnor  axes  of  ellipse  as 

shown  in  figure  1, 

»  coefficients  in  the  series  denoting  stress  for  the  plate. 

^  f  §  coefficients  in  the  series  denoting  the  stress  fonotion  for 

°      °  the  elliptic  disk. 

E    .  E  Young's  moduli  at  a  given  point  in  the  pUte  associated  nith 

^  tension  or  conQjaression  parallel  to  the  x-  and  y-azes* 

^  ,  ^        Young's  moduli  at  a  given  point  in  the  disk  associated  with 
^        tension  or  con^jression  parallel  to  the  x-  and  y-axes# 

u  ,  t'  dlspOLacements  of  a  point  in  the  plate  parallel  to  the  x- 

and  y^axes  respectively. 

S  »  displacements  of  a  point  In  the  disk  parallel  to  the  x-  and 

y-axes  respectively, 

con^onents  of  strain  in  plate  parallel  to  x-  and  yHDcss, 
coB^jonents  of  strain  in  disk  parallel  to  »-  and  y^«XBs. 
shearing  strain  for  plate  in  XT-plane, 
shearing  strain  for  disk  in  XY-^)lane, 

^  modulus  of  rigidity  associated  irith  the  shearing  strain  in 

^        the  XY-plane  of  the  plate, 

« 

yU^        modulus  of  rigidity  in  the  disk. 
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normal  components  of  stress  peirallel  to  x-  and  y-axes  in  plate. 

normal  caiiQ)onents  of  stress  parallel  to  x-  and  7-^es  in  disk* 

shearing  stress  component  for  plate  in  the  XY-plane# 

shearing  stress  co!i5)onent  in  the  Xl-plane  of  disk, 

the  Poisson's  ratio  in  the  plate  given  by  the  ratio  of  the 
contraction  parallel  to  the  y-axLs  to  the  extension  parallel 
to  the  X-axis  associated  -with  a  tension  parallel  to 
the  Xr-axis* 

the  Poisson's  ratio  in  the  elliptic  disk. 
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In  general,  a  naught  superacript  on  (or  above)  a  synbol  denotes 
relationship  to  the  elliptic  disk  iiithin  the  plate. 
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